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The longitudinal long-range correlations of the triangular flow event plane angles are calculated 
in a Boltzmann + hydrodynamics hybrid approach. The potential to disentangle different energy 
deposition scenarios is explored by utilizing two different transport approaches for the early non- 
equilibrium evolution. In the hadronic transport approach the particle production in high energy 
heavy ion reactions is mainly governed by string excitation and fragmentation processes which are 
absent in the parton cascade approach. We find that in both approaches the initial state shows a 
strong longitudinal correlation of the event plane angles which is diluted but still persists in the 
final state momentum space distributions of the produced particles. A ridge-like structure can also 
be caused by near-collinear gluon radiation in a parton cascade approach and does not necessarily 
prove longitudinal flux tubes in the initial state. 

PACS numbers: 25.75.-q,24.10.Lx,24.10.Nz 



Long- range rapidity correlations have been measured 
in heavy ion reactions at the Relativistic Heavy Ion Col- 
lider (RHIC) [2-01 ■ In two-particle correlations with a 
high pt trigger hadron in the two-dimensional At] — A.<fi 
plane a ridge-like structure around A</> = that extends 
over several units in pseudorapidity difference has been 
observed. Later on, the same phenomenon was also found 
in untriggered correlation functions and more detailed 
measurements with respect to particle composition and 
effective temperatures have been performed that indicate 
the 'bulk'-like properties of the ridge. Therefore, the orig- 
inal interpretation as a jet-medium effect has been 
questioned and other explanations as e.g. Color-Glass- 
Condensate (CGC) inspired flux tube structures that get 
boosted by radial flow have become more favored |9l-[l3j . 

Recently, initial state fluctuations have been proposed 
as the major source for many of the structures that ap- 
pear in two-particle correlations [l4l - [l7l ]. Especially the 
third Fourier coefficient of the azimuthal distribution of 
the final state hadrons in momentum space, namely tri- 
angular flow, is studied with great interest [18l - [20j . This 
new flow observable is directly related to the fluctuations 
in the initial state and is absent in hydrodynamic calcu- 
lations assuming smooth initial state profiles. Triangular 
flow is therefore independent of the collision centrality 
and very sensitive to the viscosity of the produced mat- 
ter [Hril. 

Even though a lot of progress has been made in the 
study of initial conditions, in particular regarding their 
fluctuations, it has been rather difficult to find an ob- 
servable that is directly related to the mechanism of the 
initial energy deposition. Long-range rapidity correla- 
tions seem to be perfectly suited for this purpose, since 
they need to be build up very early during the evolution 
of the reaction due to causality. In order to study the 
effect of initial state fluctuations one needs to simulate 



the whole dynamical evolution event-by-event, which is 
computationally very expensive and makes it difficult to 
study two-dimensional 2-particle correlations in detail. 

In this paper we propose a new observable that directly 
quantifies the longitudinal long-range correlation of tri- 
angular flow and explore it utilizing a state-of-the-art 
(3+l)d Boltzmann+hydrodynamics approach [25]. The 
initial conditions are either generated using hadron-string 
dynamics from the Ultra-relativistic Quantum Molecular 
Dynamics (UrQMD) approach I26l |27| or by using a par- 
ton cascade approach (PCM) [32,|33|. In both scenarios 
the longitudinal correlations of the event plane angles are 
calculated in a way that is also accessible in experiment. 

Let us first describe the two different approaches for 
the initial conditions in more detail. For the hadron- 
string dynamics, the initial bi nary nucleon-nucleon colli- 
sions are modeled in UrQMD 26-28] following the Lund 
model of nucleon-nucleon reactions [29| involving color 
flux tubes excitation and fragmentation processes that 
provide long range rapidity correlations and fluctuations 
in the energy deposition in the transverse plane. The 
other option is to decompose the incoming nuclei accord- 
ing to parton distribution functions into quarks and glu- 
ons and use a parton cascade approach to simulate the 
early non-equilibrium evolution [301 - 4331 ] . In this case, the 
cross-sections for 2— s>2 collisions and 1— >2 splittings are 
given by perturbative QCD calculations. We have chosen 
this approach, since at first glance it does not contain any 
obvious mechanism to generate long-range correlations 
and thus provides a baseline to compare to. However, as 
we shall later see, the presence of radiated parton showers 
generated after the initial hard parton-parton scatterings 
is capable of generating long range rapidity correlations. 

For Au+Au collisions at the highest RHIC energies the 
starting time for the hydrodynamic evolution has been 
chosen to be i s tart = 0.5 fm in order to fit the final state 
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FIG. 1: (Color online) Distribution of the differences of the 
coordinate space event-plane angles $3 in different rapidity 
slices in UrQMD (full line) and PCM (dotted line) initial con- 
ditions for central (b < 3.4 fm/b = 2.4 fm) Au+Au collisions 
at -Js^ = 200 GeV. 



FIG. 2: (Color online) Distribution of the differences of the 
coordinate space event-plane angles $3 in different rapidity 
slices in UrQMD (full line) and PCM (dotted line) initial 
conditions for mid-central (6 = 5 — 9 fm/b = 6.3 fm) Au+Au 
collisions at Jsnn = 200 GeV. 



pion multiplicity at midrapidity. Only the matter around 
midrapidity (\y\ < 2) is considered to be locally thermal- 
ized and takes part in the ideal hydrodynamic evolution. 
To map the point particles from the hadron/parton cas- 
cade initial state to energy, momentum and net baryon 
density distributions each particle is represented by a 
three-dimensional Gaussian distribution [34| . The width 
of the Gaussian distribution has been chosen to be a = 1 
fm for UrQMD and <r = 0.5 fm for the PCM. 

The ideal hydrodynamic evolution [Sa . [36j for the hot 
and dense stage of the collision translates the initial fluc- 
tuations in the transverse energy density to momentum 
space distributions. A hadron gas equation of state 37] 
has been used for the calculation with UrQMD initial 
conditions. 

The transition from the hydrodynamic evolution to the 
transport approach when the matter is diluted in the late 
stage is treated as a gradual transition on an approxi- 
mated constant proper time hyper-surface (see [38[ for 
details). For the hadronic calculation an energy density 
of 713 MeV/fm 3 has been chosen as a transition criterion. 
Late stage hadronic rescattering and resonance decays 
are taken into account in the hadronic cascade. 

The above event-by-event setup includes all the main 
ingredients that are necessary for the build up of trian- 
gular flow [l8[. Since the complete final state particle 
distributions are calculated, an analysis similar to those 
applied by experimentalists is used. 

In [l8j the pseudorapidity dependence of triangular 



flow in the hybrid approach explained above has been 
published. The flatness over a broad range of pseudo- 
rapidities indicates long-range correlations, but does not 
actually prove that correlation. The preferred axis of 
the triangular shape could randomly fluctuate from one 
rapidity bin to the other and one would still observe a 
flat rapidity dependence as long as the magnitude of the 
third Fourier coefficient is the same in every bin. To 
really study the long-range correlation between the 'hot 
spots' we propose to investigate the triangular flow event 
plane in different rapidity bins separately and calculate 
the correlation between these. 

In the initial state coordinate space distributions, one 
can define the event plane angles for different harmonics 
(we will concentrate onn = 2 and n — 3 in this analysis) 
in the following way: 



1 (r sm.(nd>)) 
$„ = -arctan-) ) — 

n (r n cos(n<p)) 



(1) 



where r and <f> are the polar coordinates of the particles 
in the center of mass frame of the collision. For elliptic 
flow this angle is the one that defines the so called partic- 
ipant plane. With the used conventions $„ is defined in 
the region between —ir/n and +tt /n. To explore the lon- 
gitudinal correlations, the event plane angle is calculated 
once in the whole rapidity range from — 2 < y < 2 and 
then in different bins of width |Ay| < 0.5 separately. The 
distributions of the difference between those two angles 



3 



— ' 1 1 1 1 1 1 1 - 

PCM, b=6.3 fm 



w/o rad., 1.5<y<2 

w/o rad.,1 .0<y<1 .5 

w/o rad.,0.5<y<1.0 

w/o rad.,0<y<0.5 ■ 

0<y<0.5 

0.5<y<1.0 

1.0<y<1.5 

1.5<y<2.0 




-0.5 -0.4 -0.3 -0.2 



-0.1 0.0 0.1 
$ 3 (Ay)-$ 3 



0.2 0.3 0.4 0.5 



< 

j5 3 

z 



UrQMD, b<3.4fm 
full line: a=$ 3 
dotted line: a=^ 3 



1 .5<y<2 

1.0<y<1.5 . 

0.5<y<1.0 

0<y<0.5 - 

0<y<0.5 

0.5<y<1.0 ■ 

1 .0<y<1 .5 

1.5< y<2.0- 




FIG. 3: (Color online) Distribution of the differences of the 
coordinate space event-plane angles $3 in different rapidity 
slices from the parton cascade in the default scenario (dotted 
line) and without gluon radiation (full line) initial conditions 
for mid-central (6 = 6.3 fm) Au+Au collisions at ^/snn = 200 
GeV. 



FIG. 4: (Color online) Distribution of the differences of the 
final momentum space (dotted line) vs initial coordinate space 
(full line) event-plane angles $3 in different rapidity slices in 
the hybrid approach based on UrQMD initial conditions for 
central (b < 3.4 fm) Au+Au collisions at ^/snn = 200 GeV. 



in different rapidity slices is shown in Fig. [T] We have 
checked that the results are symmetric around midra- 
pidity and the results are qualitatively unchanged if the 
bin size is varied from |Ay| < 0.25 — 1.0. Furthermore, 
the result without any additional physics mechanism that 
introduces a long-range correlation is a flat distribution 
(corresponding to a 5-function in the dN/dAy distribu- 
tion) which we can reproduce by sampling particles from 
different events together and apply the exact same cor- 
relation analysis, similar to a mixed event technique. 

In Fig. [T] and [2] the correlation of the event plane an- 
gles is shown for two different centralities and the two dif- 
ferent initial state transport approaches. In both cases, 
there is a strong correlation visible that is largest at 
midrapidity and smaller at the most forward y bin. In 
the parton cascade approach the longitudinal correlation 
can be attributed to the emission of parton showers that 
appear to be correlated to the following hard collision and 
are emitted along the beam axis, whereas in the hadronic 
transport approach the string excitation and fragmenta- 
tion mechanism offers the most plausible explanation. 

Fig. |3] shows a comparison for the parton cascade ini- 
tial conditions with and without time-like branchings, 
which initiate the parton showers. Without gluon emis- 
sion after the hard collision, the longitudinal correlation 
of the triangular flow event plane angle in the initial state 
is much smaller (if existent at all). That proves that 
the initial long-range correlation in the parton cascade is 



mainly caused by the initial gluon radiation. 

To investigate if the initial long-range correlation sur- 
vives the hydrodynamic expansion the same analysis can 
be performed for the final state momentum space event 
plane angles, 



= larctan ^;^ 
n (pt cos(n0 p )) 



(2) 



where (px, 4>p) are polar coordinates in momentum space. 
In Figs. 2] and [5] the longitudinal correlation of the final 
state event plane angles ^3 is shown in comparison to 
the previously presented distribution for $ 3 . The corre- 
lation gets significantly diluted during the ideal hydro- 
dynamic expansion, but it persists also in the final state 
particle distribution with the exception of the most for- 
ward/backward rapidity bin that is accessible in this cal- 
culation. We have checked that the same result qualita- 
tively is obtained by employing the parton cascade initial 
conditions. The final state event plane angles ^3 in dif- 
ferent rapidity bins could be measured in experiment to 
investigate the longitudinal correlations of triangular flow 
and to legitimate the assumption that the event plane 
angle does not change as a function of rapidity. 

Another way to study the initial and final state corre- 
lation of the event plane angles is shown in Fig. |6] Here, 
the differences between the two angles in different rapid- 
ity bins are shown for mid-central collisions in the hybrid 
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FIG. 5: (Color online) Distribution of the differences of the 
final momentum space (dotted line) vs coordinate space (full 
line) event-plane angles $3 in different rapidity slices in the 
hybrid approach based on UrQMD initial conditions for mid- 
central (b — 5 — 9 fm) Au+Au collisions at ^/snn = 200 GeV. 



approach based on UrQMD initial conditions. As a com- 
parison the correlation for elliptic flow at is also shown. 
For the three rapidity bins between and 1.5 the correla- 
tion is clearly visible for triangular flow as well, whereas 
it is less pronounced in the most forward/backward bin. 
This confirms the conclusion of the Fig. [5] that the long- 
range correlation only survives over three units in rapid- 
ity around y = 0. 

In summary, we have proposed a new way to inves- 
tigate the long-range correlations of triangular flow in 
heavy ion collisions. A hadron and a parton cascade ap- 
proach both lead to sizable long-range correlation in the 
initial state. These correlations are translated to the final 
state hadron distributions, but are significantly diluted 
and smoothed out during the final hadronic rescattering. 
Measuring the triangular flow event plane angles in dif- 
ferent rapidity bins serves as a pre-requisite to legitimate 
the assumption of a constant event plane over rapidity 
as it is commonly assumed by experimental groups. The 



second conclusion is that one needs to take into account a 
realistic dynamical evolution (e.g. [39|,[4(|) to prove that 
the long-range correlations that might be initially estab- 
lished also survive to the final state particle distributions 
instead of simple parametrizations. 
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FIG. 6: (Color online) Event-by-event correlation between 
initial and final state event plane angles for elliptic (dotted 
line) and triangular flow (full line) in different rapidity slices 
in the hybrid approach based on UrQMD initial conditions for 
mid-central (b = 5 — 9 fm) Au+Au collisions at ^/snn = 200 
GeV. 
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